Chapter 1 

Understanding Hawking radiation from 
simple models of atomic Bose-Einstein 
condensates 
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Abstract This chapter is an introduction to the Bogoliubov theory of di- 
lute Bose condensates as applied to the study of the spontaneous emission 
of phonons in a stationary condensate flowing at supersonic speeds. This 
emission process is a condensed-matter analog of Hawking radiation from 
^ I astrophysical black holes but is derived here from a microscopic quantum 

theory of the condensate without any use of the analogy with gravitational 
^ [ systems. To facilitate physical understanding of the basic concepts, a sim- 

bX) pie one-dimensional geometry with a stepwise homogenous flow is considered 

which allows for a fully analytical treatment. 



1.1 Introduction 

One of the most spectacular predictions of Einstein's General Relativity is the 
existence of Black Holes (BHs), mysterious objects whose gravitational field is 
so strong that not even light can escape from them but remains trapped inside 
a horizon. According to the standard view, BHs arc formed by the collapse 
of massive stars (M > 3Msun) at the end of their thermonuclear evolution 
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when the internal pressure is no longer able to balanee the gravitational self 
attraction of the star. Furthermore supermassive BHs {M > lOMgun) are 
supposed to constitute the inner core of active galaxies. 

As no light can escape from them, BHs are expected to be really "black" 
objects. In particular, their observational evidence can only be indirect: typ- 
ically, the presence of a black hole is deduced by observing the behavior of 
matter (typically hot gas) orbiting outside the horizon. A hypothetical iso- 
lated BH (i.e. a BH immersed into vacuum) would not manifest its presence 
except for his gravitational field, which after a short transient time becomes 
stationary (even static if there is no angular momentum) . 

In 1974 Hawking showed [1] that this common belief is incorrect. If one 
takes into account Quantum Mechanics, static and stationary BHs are no 
longer "black" , but rather emit a steady radiation flux with a thermal spec- 
trum at a temperature given, simply speaking, by the gradient of the gravita- 
tional potential at the horizon. This intrinsically quantum mechanical process 
is triggered by the formation of the horizon and proceeds via the conversion 
of vacuum fluctuations into on shell particles. This effect is a universal feature 
of BHs, completely independent of the details of the BH formation. 

In spite of the interest that this fascinating effect has raised in a wide 
audience, no experimental evidence is yet available in support of this amazing 
theoretical prediction. Since the emission temperature scales as the inverse 
of the BH mass (T - 10"'^ K for a solar mass BH), the expected Hawking 
signal is in fact many order of magnitudes below the 2.7 K of the cosmic 
microwave background. As a result, the Hawking radiation by BHs appears 
to be a completely irrelevant process in any realistic astrophysical situation, 
with no hope to be detected in the sky. This situation is rather frustrating, 
since the conceptual relevance of Hawking discovery is extremely profound: 
the existence of Hawking radiation Hawking radiation allows such a beautiful 
synthesis between gravity and thermodynamics that it cannot be just an 
accident; many people indeed regard Hawking result as a milestone of the 
still to be discovered quantum theory of gravity. 

After almost 40 years of research on BHs, the attitude nowadays appears 
a bit different and more promising on the experimental side. In particular, it 
was realized that the Hawking emission process is not at all peculiar to grav- 
itational physics: its "kinematical" rather than "dynamical" nature makes 
it manifest itself in different physical contexts. This way of looking at the 
Hawking effect has its origin in a paper by Unruh in 1981 [2j where a steady 
emission of thermal phonons was predicted to appear in any fluid in stationary 
flow turning supersonic: the basic process underlying this phonon emission 
is completely identical to the one discussed by Hawking for the gravitational 
BH, in the sense that the mathematical equations describing it are exactly 
the same as the ones describing Hawking radiation from gravitational BHs. 
The reason for this amazing and unexpected "analogy" is that the equation 
describing the propagation of long wavelength sound waves in a moving fluid 
can be recast in terms of a massless scalar field propagating in a curved space- 
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time with a suitably chosen "acoustic metric" . In particular, the point where 
a sub-sonic flow turns supersonic plays the role of an "acoustic horizon" , 
since sound waves in the supersonic region are no longer able to propagate 
upstream. As it happens to light inside a BH, sound waves are trapped inside 
the sonic horizon of the "acoustic black hole" : upon quantization, it is then 
straightforward to expect the emission of Hawking radiation by the horizon. 
Nowadays, we know that this analogy with gravitational systems is not lim- 
ited to fluids but can be developed for many other condensed matter and 
optical systems Unlike gravitational BHs, these condensed matter analog 
models often possess a well understood quantum description at the micro- 
scopic level, which allows for a complete control of their physics. This is the 
case, in particular, of atomic Bose-Einstein condensates which are the subject 
of the present chapter. 

The relevance of the analogy is therefore twofold. At one hand, one can 
concretely consider investigating the actual existence of Hawking radiation 
using table top experiments with a complete control of the physical system. 
On the other hand, the detailed knowledge of the underlying microscopic 
quantum theory of these systems allows us to address a very delicate point 
in the theory of Hawking radiation and possibly to eliminate some intrinsic 
inconsistencies of its standard derivation. 

In the absence of a complete and self consistent quantum theory of grav- 
ity, one typically adopts a semiclassical framework where gravity is treated 
classically according to General Relativity, whereas light and matter fields 
propagating on top of the curved space time are quantized. This is the so 
called Quantum Field Theory in Curved Space [3]. One expects this scheme 
to provide a sufficiently accurate description of the gravity-matter systems 
for scales sufficiently large as compared to the fundamental quantum scale 
for gravity, the so-called Planck scale equal to 10"^'^ cm or 10^^ GcV. Ap- 
proaching this Planck scale, one can reasonably expect that this semiclassi- 
cal description becomes inaccurate and has to be replaced by a (yet to be 
discovered) complete theory of quantum gravity. 

Now because of the infinite (exponential) redshift suffered by the Hawking 
phonons in their journey from the near horizon region to infinity, a given 
mode of Hawking radiation measured at time t with frequency v far from the 
BH appears to have had a frequency v' = v e'^*^^^ near the BH horizon (i? is 
the radius), which rapidly overhangs Planck energy: this feature makes the 
derivation clearly inconsistent and casts serious doubts on the very existence 
of Hawking BH radiation. This is the so called transplanckian problem [5] . 

The same kind of argument can be repeated also for Hawking like radi- 
ation in condensed matter systems: because of the infinite Doppler shift at 
the sonic horizon, the modes responsible for Hawking like radiation oscillate 
near the horizon at a wavelength much smaller than the intermolecular or 
interatomic spacing, which makes the hydrodynamical long wavelength ap- 
proximation inconsistent. On this basis, it would therefore be difficult to rule 
out the possibility that Hawking radiation be an artifact due to an illegiti- 
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mate extrapolation of the long wavelength approximation, that is a spurious 
outcome without any physical reality. 

In this perspective, analogue condensed matter systems provide a new 
angle from which the transplanckian problem of Hawking radiation can be 
attacked: as they possess a detailed and well understood microscopic quan- 
tum description, the question of the existence of Hawking radiation can be 
investigated from first principles, without any use of the hydrodynamical ap- 
proximation and hence of concepts borrowed from the gravitational analogy 
like effective metric, horizon, etc. So far, most of the work in this direction 
has been performed using atomic BECs, but the fully positive answer coming 
from these studies appears to hold under very general assumptions: Hawking 
radiation is indeed a real physical phenomenon! 

A closer look at the spectral and coherence properties of the predicted 
Hawking radiation match the original expectation that, if the transition is 
sufficiently smooth with respect to the microscopic scales of the fluid, the 
Hawking emission of Bogoliubov phonons is thermal at a temperature pro- 
portional to the gradient of the flow potential at the sonic point [6]. In ad- 
dition, several novel interesting features have pointed out in regimes beyond 
the hydrodynamical approximation as well as in different conflgurations, e.g. 
white holes (the time-reversed of a black hole) [7] and the so-called black-hole 
lasers (a pair of adjacent black and white hole horizons) [27] . 

Simultaneously to these theoretical and conceptual advances, a great effort 
has been devoted in the last years to the identiflcation of the most promising 
physical systems where to experimentally investigate the analogue Hawking 
radiation. 

Having established that the Hawking effect exists, one can start to think at 
the best experimental setting to reveal it. There are many systems proposed 
at this end, like ultracold atoms, optical systems, water tank experiments and 
others. At the moment, experiments with water tanks [8] have detected the 
classical counterpart of Hawking emission in flows showing white hole hori- 
zons: stimulated emission by the Hawking mechanism is probed by sending 
a classical incident wavepacket of surface waves against the horizon. Unfor- 
tunately, these experiments at room temperature do not appear suitable to 
investigate the very quantum phenomenon of Hawking radiation, that is the 
conversion of zero-point fluctuations into observable quanta by the horizon. 
An observation of Hawking radiation from laser pulses propagating in nonlin- 
ear optical media has been recently reported [9] , but this result is still object 
of intense discussion in the community |10j . 

The main experimental difficulty in the quest for analog Hawking radiation 
in condensed matter systems is the extremely weak intensity of the signal in 
realistic systems, which makes it to be easily covered by competing effects 
like the thermal emission due to the non zero temperature of the systems as 
well as quantum noises. In this respect, atomic gases appear as most promis- 
ing systems |13] . as they combine a variety of tools for the manipulation of 
the diagnostic of the system down to microscopic level, to the possibility 
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of cooling the system to very low temperatures where the zero point quan- 
tum fluctuations start playing an important role. Still, even in these systems 
temperatures lower than the expected Hawking temperature of the order of 
10 nK are hardly reached, and further difficulty comes from the detection of 
the Hawking phonons emitted from the horizon. 

A major breakthrough that appears to bypass both these problems was 
proposed by us in 2008 [11] and is based on the use of density correlations , a 
modern powerful tool to investigate microscopic properties of strongly corre- 
lated atomic gases and in particular of their elementary excitations. Taking 
advantage of the fact that the Hawking radiation consists of correlated pairs 
of quanta emitted in opposite directions from the horizon, a characteristic 
signal will appear in the density-density correlation function for points sit- 
uated on opposite sides with respect to the horizon. This unique signature 
was made quantitative using gravitational physics methods and then numer- 
ically confirmed by ah initio simulations of the condensate dynamics based 
on a microscopic description of their collective properties [12] . As a result, it 
appears to be an ideal tool to isolate the Hawking radiation signal from the 
background of competing processes and of experimental noise even at non- 
zero temperatures. Of course, a similar strategy would be clearly impossible 
in astrophysical black holes, as no access is possible to the region beyond the 
horizon. 

In this paper we shall use standard tools of the theory of dilute Bose 
gas to show in a rather pedagogical way how Hawking radiation emerges 
in an atomic BEG and to explain its features on a simple and analytically 
tractable toy model. Our treatment, as we shall see, resembles very much 
learning elementary Quantum Mechanics using one dimensional Schrodinger 
equation with square potentials. Most of the material presented here was 
originally published in Refs.[18] and [T7] . 



1.2 The theory of dilute Bose-Einstein condensates in a 
nutsheU 

In this section we give a brief and rapid introduction to the theory of BECs. 
In particular, we shall review the Gross-Pitaevskii equation describing the dy- 
namics of the condensate at mean field level and the Bogoliubov description of 
quantum fluctuations of top of it. More details can be found in textbooks [25] 
and in dedicated reviews [29] . 

Bose-Einstein condensation is characterized by the accumulation of a 
macroscopic fraction of the particles into a single quantum state. To achieve 
such a quantum degeneracy very low temperatures are required (on the order 
of T = 100 nK for the typical densities of ultracold atomic gases in magnetic 
or optical traps) , where particles are no longer distinguishable and their Bose 
statistics start being relevant. 
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1.2.1 The Gross-Pitaevskii equation and the 
Bogoliubov theory 

The model Hamiltonian describing a many-body system composed of N in- 
teracting bosons confined in an external potential Vgxti^) can be written in 
a second quantized formalism as: 



H = d^x 



(1.1) 



where !J'(t, x) is the field operator which annihilates an atom at position x 
and obeys standard bosonic equal time commutation rules 

[<f-(x),if-^(x')] = (53(x-x'). (1.2) 

The model Hamiltonian (|l.ip is generally used within the dilute gas approx- 
imation where the two body interatomic potential can be approximated by 
a local term V{x — x') = ^^^(x — x') with an effective coupling constant g 
related to the atom-atom scattering length a by g = AiTli^a/m. 

At sufficiently low temperatures well below the Bose-Einstein condensation 
temperature, a macroscopic fraction of atoms is accumulated into the single 
one-particle state of lowest energy, described by the macroscopic wavefunction 
!fb(x). The time evolution of the macroscopic wavefunction in response of 
some excitation (e.g. a temporal variation of the confining potential Vext) is 
described by the Gross-Pitaevski equation 

whose form can be heuristically derived by performing a mean-field approxi- 
mation 'F ^ in the Heisenberg equation 



^(i,x),i? (1.4) 



for the time-evolution of the atomic quantum field operator F. The ground 
state wavefunction naturally emerges as the lowest-energy steady-state ^'o(x) 
of the Gross-Pitaevskii equation and oscillates at a frequency fj,/h. 

Small fluctuations around the mean-field can be studied within the so- 
called Bogoliubov approximation, where the bosonic field operator ^ is writ- 
ten as the sum of a mean-field classical field plus quantum fiuctuations. In its 
usual formulation to describe weakly excited condensates, one takes a steady 
state !fo as the mean-field, 

#(<, x) - tf'o(x) [1 + kt, x)] e-'"*/^ (1.5) 
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The field operator 4> describing fluctuations then satisfies the Bogoliubov-de 
Gennes (BdG) equation 

in-^ = -(^— + ^vU + ng + , 1.6 

at \ 2m m Wq J 

where n = |!^bP- The next subsections will be devoted to a rewriting of 
the BdG equation in terms of a curved space-time with an effective metric 
determined by the spatial profiles of the local speed of sound c = y^ng/m 
and of the local flow velocity Vq. 



1.2.2 Analogue gravity in atomic BECs 

We stop for the moment the formal development of BEG theory and show 
how a different parametrization of the field operator leads to a reinterpreta- 
tion of the above equations in a hydrodynamical language and then to the 
gravitational analogy [3] . 

Using the so called density-phase representation of the condensate wave- 
function if'o = \fne^^ ^ the Gross-Pitaevskii equation (|1.3p can be rewritten 
as a pair of real equations, 



9*" + V(nv) =:0 , (1.7) 
2m" 



h.dtQ = - — {'^6?-9n-Ve.t-y,i- (1-8) 



the former equation Eq. (|1.7|) is the continuity equation with an irrotational0 
condensate velocity vq = hVO/m. The latter is analogous to Euler equation 
for an irrotational inviscid fluid, with an additional "quantum pressure" term 

K,(x) 

Zm \Jn 

describing a kind of stiffness of the macroscopic wavefunction. 

In this density-phase representation, the Bogoliubov expression (|1.5p of 
the fleld operator is rewritten as 



and the Bogoliubov equation (jl.6|) reduce to a pair of equations of motion 
for the fluctuations in the density h\ and in the phase (Q\) in the form 



^ From the definition of the velocity field vq, it is immediate to see that the vorticity 
in the condensate can only appear at points where the density vanishes. 
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hdA = -hwoV§i - — ni + ^^^■'vlnvC-^)] = , (1.11) 
fill 

dthi = -V(vo7ii + —Vei). (1.12) 

m 

Here, a fundamental length scale is set by the so-called healing length defined 
as ^ = h/mc m terms of the local speed of sound c = ^Jngjrn. 

If one is probing the system on length scales much larger than ^ (the 
so-called hydrodynamic approximation ), the last term in eq. (jl.lip can be 
neglected. As a result, the density fluctuations can be decoupled as 



hn 

ni = TT 



(1.13) 



When this form is inserted in eq. (|1.12p , the equation of motion for the phase 
perturbation 

n Tt 
- (dt + Vvo) — ^{dt + voV)0i + V-V01 = (1.14) 

can be rewritten in a matrix form 

d^ir^dj^) = (1.15) 
where the matrix elements /^"^ are defined as 

r = -4, = = -^vl, P = ^^{cH'^ - vlvi) (1.16) 

in terms of the condensate density n and local velocity vg. Greek indices 
jji^v ~ 0,1,2,3 indicate 4-dimcnsional objects, while latin ones i — 1,2,3 
indicate the space coordinates. 

Now in any Lorcntzian manifold the curved space scalar d'Alembertian 
operator can be written as 

□ = -^d^i^g^'^d,) (1.17) 



where g is the metric, g^^ its inverse and g = det{g^y). Keeping this in mind, 

Eq. (|1.14[) for the condensate phase dynamics can be rewritten in the form 
of a curved space wave equation 

□6*1=0, (1.18) 

provided one identifies 

^gg^^'^r, (1-19) 
which can be inverted leading to the effective metric 
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In summary, we have shown that under the hydrodynamical approxima- 
tion, the equation of motion for the phase fluctuation in a BEC can be rewrit- 
ten in terms of a Klein-Gordon equation for a massless scalar field propagating 
in a fictitious space-time described by the metric 17^1/ defined by Ea. (|1.20p . 
This is the core of the gravitational analogy . 

One should stress that this Lorentzian space-time has nothing to do with 
the real space-time in which our BEC lives. Note also that the invariance of 
eq. (|1.18p under general coordinate transformation is fake. The underlying 
BEC theory is not even (special) relativistic, but Newtonian, with an absolute 
time, the laboratory time, with respect to which the equal time commutators 
(eq. (|1.2p ') are given. 

Said this, one can give a closer look at the metric g^i, given by Ea. (|1.20p : a 
particularly interesting situation is when a steady fluid turns supersonic (i.e. 
I Vol > c) in some region of space. In a gravitational analogy , such a con- 
figuration corresponds to a black hole as described in the so-called Painleve 
-Gullstrand coordinate system and is therefore called a "sonic black hole" : as 
sound waves travel at a velocity c lower than the fluid velocity vq, they are 
not able to propagate back and result trapped inside the supersonic region 
beyond the "sonic horizon", i.e. the locus where |vo| = c. 

In such a setting Eq. (ll.lSp describes a massless scalar field propagating in 
a black hole space-time. But this is exactly the system considered by Hawk- 
ing to obtain his famous result. One can then repeat step by step Hawking's 
derivation of black hole radiation. First of all, one has to expand the field in 
modes and focus his attention on those upstream propagating modes which 
are barely able to avoid being trapped by the horizon formation and escape 
in the subsonic region. Upon quantization, comparison of the 'in' and 'out' 
vacuum states then shows that they are inequivalent since the corresponding 
annihilation and creation operators are related by a Bogoliubov transforma- 
tion that mixes them in a non-trivial way. As a result, one can expect that 
an emission of Bogoliubov phonons by the horizon appears in the sub-sonic 
region, which are thermally distributed at a temperature given by the surface 
gravity k, of the sonic horizon defined as 



1 d{c 

with n — 



2c dn 



(1.21) 

hor 



where n is the spatial coordinate normal to the horizon. 

It is however crucial to keep in mind that this conclusion is based on a 
very strong assumption, namely the long wavelength approximation, which 
has allowed to neglect the last term in eq. (|l.lip . to rewrite this equation 
as D^^i = 0, and to introduce the gravitational analogy. As explained in 
the introduction, the modes of the field responsible for the Hawking emission 
experience an infinite Doppler shift when leaving the near horizon region in 
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the upstream direction. As their wavelength in this region is many order of 
magnitude smaUer than the heahng length of the atomic gas, all the derivation 
of Hawking radiation in atomic BEC outlined above is at least questionable. 

For this reason we go back to the original microscopic BEC theory of 
Sec ll.2l and try to derive Hawking radiation without making any hydrody- 
namical approximation and without any reference to the gravitational analog: 
the emission of Hawking radiation in BEC supersonic configurations will then 
appear as a natural outcome of the underlying quantum theory. 



1.3 Stepwise homogeneous condensates 

A simple analytical treatment can be developed to show the occurrence of 
Bogoliubov phonons creation " a la Hawking" in an atomic BEC undergoing 
supersonic motion in a very idealized setting consisting of two semi-infinite 
stationary homogeneous one dimensional condensates (left and right sector) 
connected by a step- like discontinuity [TH|. As the purpose of this article is 
mostly a pedagogical one, we do not enter into a discussion of the actual 
experimental feasibility of this configuration, for which we refer to the most 
recent research literature [30l [31] . 

In particular, we assume the condensate to have a everywhere uniform 
density n in both sections as well as a spatially uniform flow velocity v along 
the negative x axis. The external potential Vext and the repulsive atom-atom 
interaction coupling g are supposed to be constant within every sector, but 
to have different values in each sector, satisfying 

K',t+.g'n = K';t+/n. (1.22) 

Here, the superscripts "1" and "r" refer to left {x < 0) sector and right {x > 0) 
sector respectively, the discontinuity being located at a; = 0. The change in 
the interaction constant g can be obtained cither via the dependence of the 
atom-atom scattering length on the value of a static external magnetic field, 
or by modulating the transverse confinement orthogonal to the x direction. 
Such a change in g directly reflects onto the local sound speed, that then has 
different values c' and in the two sections, defined as usual by m(c'''')^ = 
ng'"'^ . Thanks to the condition (|1.22p . the plane wave form 

tZ'o(i,x) V^e*'^'°^-''"°* (1.23) 

of the condensate wavefunction is a solution of the Gross-Pitaevski equation 
(|1.3p at all times t and positions x. The wavevector kg and the frequency ujq 
are related to the flow velocity vq by vq = hko/m and hujo = fi^k^/ (2m) +gn. 

Let us look now at the solutions of the BdG equation Eq. (|1.6p for the 
fluctuation fleld within each sector. Exploiting the stationarity of the con- 
figuration, it is convenient to split (j) into its "particle" and "antiparticle" 
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(1.24) 



where aj and aj are phonons annihilation and creation operators, satisfying 
the usual bosonic commutations rules [di, dj] = 6ij . The mode functions 
and ifj satisfy the motion equations 



(1.25) 



that follow from Eq. (|1.6p and its conjugate and are to be chosen as oscillating 
at a frequency LOj . Imposing that the equal time commutators satisfy 



[(i>{t,x),(t>''{t,x')] = -6ix-x') 



(1.26) 



provides the normalization of the modes 



hn 



(1.27) 



The sum over j in (|1.24[) only involves positive norm modes for which the 
sign in p.27p is positive. 

Within each of the two a; < and x > spatially uniform regions, the 
mode functions have a plane wave form 



(1.28) 



where D{uj) and E{uj) arc normalization factors to be determined using cq. 
(fL?f)) . Inserting eqs. into yields 



{ll> - vok) - 
(uj - Vok) 



^ck'^ c 



(1.29) 



the existence of nontrivial solutions requires that the determinant associated 
to the above homogeneous system vanishes. 



{uj - voky =c^\k^ + 



(1.30) 



Solving this implicit equation provides the so-called Bogoliubov dispersion of 
weak excitations on top of a spatially uniform condensate, 
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Vok = icy 



f2±{k) : 



(1.31) 



here, f2± is the excitation frequency as measured in the frame co-moving with 
the fluid. The + (— ) sign refers to the positive (negative) norm branch. As 
expected, for small k such that fc^ ^ 1, the dispersion relation is linear 



VQk = ±ck 



(1.32) 



this is the hydrodynamical regime to which the gravitational analogy is 
strictly speaking restricted. At higher k, the corrections to the linear disper- 
sion are positive and the modes propagate supersonically. For large k such 
that k^ ^ I, the relation tends to the quadratic dispersion of single particles. 
The normalization condition gives 



\D{cu)\'-\E{co)\'^± 



1 



2nhn 



dk 



duj 



(1.33) 



which using eqs. (|1.29p yields the normalization factors 

D{u;) = 



vok + ^ 



ATThnc^k"^ 



E{u) 



(---ofc)(£: 



^A^hncik^\{uj~vok){^) ' 



(1.34) 



as expected, positive (negative) norm states correspond to the branch of the 
dispersion relation at a positive (negative) comoving frequency. Remarkably, 
for any positive norm branch of frequency tu and wavevector k, there exists a 
negative norm branch of opposite frequency —oj and wavevector —k. Taking 
advantage of this duality, one can use both positive and negative norm states, 
replacing the sum over j in (jl.24p with an integral over cj and restrict to the 
positive frequency ones. 

Let us go back to the dispersion relation eq. (jl.SOp . At fixed w (> 0) this 

(i) 

is a fourth order equation in k. It admits four solutions fc^I) and in general 
(jji^ is a linear combination of four plane waves constructed with the k's of 
the form 



(1.35) 



where the Ai(uj) are the amplitudes of the modes, not to be confused with 
the normalization coefficients D{oj). Similarly for ip^{x,t). 

As said before, our systems consist of two semi-infinite homogeneous con- 
densates joined at x=0 where there is a step-like discontinuity in the speed 
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of sound. Looking at the modes equations eqs (jl.25p , one has to require that 
the solutions in the left region and the ones in the right region satisfy the 
following matching conditions at x = 

[0] = 0, [0'] = 0, = 0, y] = 0, (1.36) 

where [f{x)] = lime_j,o[/(a; + e) ^ f{x — e)] and ' means These four condi- 
tions allow to establish a linear relation between the left and right amplitudes 

A\ = A-U,A) (1.37) 

where M is a 4 x 4 matrix called the matching matrix, not to be confused 
with the scattering matrix S we will introduce later, whose dimensionality 
may vary. 

To proceed further in the analysis and explicitly solve the dispersion re- 
lation to get the four roots /ci*-*, one has to specify the flow configuration 
under investigation, as the position of the roots in the complex plane varies 
according to the subsonic or supersonic character of the flow. In the next 
sections we shall separately consider the different cases. 



1.4 Subsonic-subsonic configuration 

We start by considering a flow which is everywhere sub-sonic , that is with a 
flow speed vq smaller in magnitude than both q and Cr, that is |wo| < c^.;. A 
sketch of the configuration under investigation is given in the upper panel of 
FigHH 



The Bogoliubov modes and the matching matrix 

The Bogoliubov dispersion in a subsonic flow is graphically displayed in the 
two lower panels of Fig ll.ll for two different relative values of the sound speed 
c' (left) and d' > c' (right). The positive (negative) norm branches are plotted 
as solid (dashed) lines. For any given a; > 0, two real solutions belonging to 
the positive norm branch exist within each I, r region: one, ku, has a positive 
group velocity Vg = ^ and propagates in the rightward, upstream direction; 
the other, fc„ has a negative group velocity and propagates in the leftward, 
downstream direction. The u, v labels used to indicate these solutions are 
the conventional ones in General Relativity. These two real roots admit a 
pcrturbative expansion 
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Fig. 1.1 Upper panel: sketch of the subsonic-subsonic flow configuration. Low panels: 
dispersion relation of Bogoliubov modes in the asymptotic flat regions away from the 
transition region. 



vq - c \ 8[vq - c)-^ J 
ku = ^ U "'Y ,3 + Oiz^)) (1.38) 

where the dimensionless expansion parameter is 2 = /c. To zeroth order in 
z, one recovers the well known hydrodynamical results fc„ = w/(wo — c) and 
fcu = w/ (wo + c). In the following, we shall indicate as ku,v the value of these 
roots in each of the two homogeneous sections on either side of the interface. 

The other two solutions of the dispersion relation are a pair of complex 
conjugate roots. Within the right sector at a; > 0, we call /cl the root with 
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positive imaginary part, which represents a decaying mode when one goes 
away from the horizon in the positive x direction. The other solution with a 
negative imaginary part corresponds instead to a growing (and therefore 
non-normalizable) mode. The opposite holds in the left sector at a; < 0; the 
k\_ root with a positive imaginary part represents a growing mode away from 
the horizon, while the other root fcL with a negative imaginary part represents 
the decaying mode. Within each Z,r region, the wavevector of these modes 
can be expanded in powers of z = C^^/c as 



0^ 



4(c2 - ^2)3 



± 



(c^ + 2vl)c^z'' 

8(c2 - 1,2)3 



(1.39) 



In summary the decomposition of (j)^ and (p^j in the left (right) regions 
reads 



l(T) 



A'1''D 



0-AO) 



We stress again the fact that the coefficients A''^'^^ j_ are the amplitudes of the 

different modes, not to be confused with the normalization coefficients, D^^^^ ^ 

for and E^Jj_ for ipi^: these latter are uniquely fixed by the commutator 
relations and the equation of motion, while the amplitudes depend on the 
choice of basis for the scattering states as we shall see in Sect. 11.4.21 Note 

and ipi^ as required by the 



that the amplitudes A^^^^j. are the same for 
equation of motion. 

The matching conditions at a: = 0, [(p] = 0, [ip] = 0, [(/)'] = 0, 



= 



impose a linear relation between the four left amplitudes A[^ ^ j_ and the right 



ones A;_„,± 

\ALJ 

where the 4x4 matrices Wi(^r) are 



Wi 



Wr 



( Al\ 

Al 

A+ 
\A-J 



(1.42) 
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/(r) 



W, 



l{r) 



D 



l{r) 



( dT^ dT^ d 



l{r) 



E 



Multiplying both sides by ^ one finally gets 



(1.43) 



A-ii 






= M 


K 






\A-I 



(1.44) 



in terms of the matching matrix M = Wi ^Wr whose explieit form is rather 
involved and is not given here. 



1.4.2 The "in" and "out" basis 





Fig. 1.2 Sketch of the Bogoliubov modes involved in the "in" (left panels) and "out" 
(right panels) basis. The mode labels refer to the dispersion shown in the lower panels 
of FigO 
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Wc now proceed to construct a complete and orthonormal (with respect 
the scalar product eq. p.27|) ') basis for the scattering states of our operator. 
This can be done in two ways: either choosing a "in" basis constructed with 
incoming scattering states (i.e. states that propagate from the asymptotic 
regions x = ±00 towards the discontinuity at a; = 0) or an "out" basis 
constructed with outgoing scattering states (i.e. states that propagate away 
from the discontinuity to x — ±00). 

Let us start with the "in" basis, whose construction is sketched in Fig jl.2l 
We define the in ii-mode (/>^'™ as a scattering state with a unit initial am- 
plitude on the left moving v mode incident on the discontinuity from the 
right (x = +00), i.e. D^, 5"'"*+''"'"'^. The incident wave is scattered by the 
discontinuity at a; = into a transmitted u-mode in the left region with am- 
plitude A^^ (i.e. A\D\ g-iut+ik^x--^ g^^^ partially reflected in the right region 
with amplitude A^ (i.e. A^^D^ ^-u^t+ik^xy order to complete the construc- 
tion, we have to include in both regions the complex decaying modes as well: 
A!j_Z)!i^e~*'^*+*'^+^ and ^'_£)'_e-*'^*+''=-^. Growing modes are not included as 
they diverge at infinity. 

The general matching equation (jl.44p becomes in this case 












= M 


K 





K 


\Al) 




\ 1 



(1.45) 



Treating M pcrturbatively in zi = we obtain (for the simplest case fo = 0; 
for the general subsonic f 7^ case the amplitudes are given in the appendix 

of HZ]) 



-A,, 



T 



2y/QC; {ci - Cr) "^Zi Cl [ci - Cr) ^ (cf C^) z} 



+ Cr cV^ (ci + Cr) 



2c3 (q ^CrY 



(1.46) 



Al 



ici {Cl - Cr) ^Zl Cl {Ci - Cr) {2cf - 3cfCr + 2c;C^ + C^) zj 



Cl - Cr 
Cl + Cr cl {Cl + Cr) 

{Cl - Cr) -Jzi {ci 



(q + Cr) 2L'Lcf/^ (q + Cr) 
Cl {-Cl + Cr) ^l ^ cf {Cl - Cr) zf 

2DlcV^ (q + Cr) 



4cf, {Ci + Cr) 
[Cr"^ + i (C(^ + Cr^ 



HA7) 



i?;c?/' {Ci+Cr) 



CrCl)] , (1.48) 

[ci+i{ci-2cr)] . (1.49) 



Note that these combine in such a way that the unitarity relation |i?p + |rp ~ 
1 is satisfied. Note also that even if they do not enter the unitarity relation, the 
amplitudes of the decaying modes are part of the full mode and their presence 
will show up explicitly contributing to the density-density correlation as we 
shall see. 
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In a similar way we can construct the (^^'*" as a scattering state with a 
unit initial amplitude on the right moving u mode incident on the disconti- 
nuity from the left (x = — oo), which is partially reflected back and partially 
transmitted, as shown in Fig ll.2l Here too, we have to include the decaying 
modes. 

The matching relation now reads 



(1.50) 







( ' \ 


1 















\ / 



yielding 

Al^R' 

Al 



Cr - Cl i {Cl - Cr) Zl ^ (q - Cr) (cf + 2cfCr 
Cl + Cr 



Cr {ci + Cr) 
.^l{ci - Cr) '^Zi 



4c3 {Cl + Cr) 

ci {Cl - Cr) 2 (c^ - 4QCr + (?r) Z, 



-(1-51) 

2 
I 



Cl + Cr 
{Cl - Cr) 



Cr [Cl + Cr ) 



ScV^ {Cl 



Cr) 



(ci - Cr) 



Diy/ci (ci + Cr) 2DLy/ciCr {Cl + Cr] 
^/ci {-Cl + Cr) y/zi y/ci{ci-Cr) 
DWr{ci+Cr) 2DWf (q + Cr) 



+ i (2q - Cr)] zf 



(1.52) 
(1.53) 



2 ■ { 1 2 
Cl + I (Ci + Cr" 



.)](i-?54) 



which implies + = 1, as required by unitarity. 

The scattering modes </>^'™ and (/)^'*", and the similarly constructed <p^'*" 
and </?Ji'™, constitute a complete "in" basis for our field operator 0, that can 
be then expanded as 

/•oo 

<^(x,t)= / a^'™<^™(t,x)+a^'"C;(^,^)+«::'"V™*(i,^)+«"'"^ 
Jo '- 

(1.55) 

The "in" vacuum |0, m) is defined as usual by aj^'*"|0, m) = and a^'™|0, in) = 
0. The A^-phonons states that constitute the "in" basis of the Hilbert space 
are constructed by a repeated action of the creation operators aj^"'™ and 
-tu.jn |-]-^g vacuum state. 

While the "in" basis has been constructed using incoming scattering 
modes, an alternative "out" basis can be constructed starting from the out- 
going scattering basis of the field operator, composed of modes that emerge 
from the scattering region around a; = with unit amplitude on a wave prop- 
agating at t = -l-oo either rightwards towards x = +00 or leftwards towards 
X = —00. 

We begin by defining the </>Ji'°"* scattering mode: as it is sketched in 
Fig. II. 2[ this is a linear combination of in-going right and left moving compo- 
nents with amplitudes A[^ and A!^ and decaying modes with amplitudes AL 



1 Simple models of Hawking radiation in atomic BECs 



19 



and A^. These coefficients are chosen in a way to give after scattering only 
a left moving w-mode of unit amplitude. This imposes the condition: 



(1.56) 



( \) 








= M 










\al) 




\ J 



that yields 

Al = R' 



Ci + Cr Cr {Cl + Cr) 



{Cl - Cr) ^Zi ^ (q - Cr) {cf + 2cf - 3qC^ + '^^'^fk^ ^^ 



4cf {Cl + Cr) 



r _ rr.,* _ , (ci - Cr) ^Zi {Cl - Cr) ^ {cf - AciCr + cf) zf 

A„ = 1 — 1 — Yjj- U-3oj 

8Cr (Cl+Cr) 

[Cr + I i2ci ~ Cr)] , (1.59) 



A' 



Cl + Cr cV {Cl + C,.) 

(Ci-Cr)^l (Cl-Cr). 



D^y/Fl {Cl + Cr) 2DLcr{ci+Cr) 



Q {-Cl + Cr) yJTl (ci - Cr) zf 



D\Cr [ci + Cr) 2Dlcl {ci + Cr) 



I (Cl + Cr 



- c.(4)<^0) 



with + |r'*|2 = 1. 

The same procedure can be used to construct the mode 0"'°"*, by imposing 
the out-going waves to consist of a unit amplitude right moving u-modc only. 
In this case, the matching relations are 





\ALJ 



M 



1 

a; 
V / 



(1.61) 



with 

Ai. = T* 



2^fc{c; iy/ci{ci - Cr)'^Zi ^l {ci - Cr) ^ (cf - 4ciCr + C^) zf 



8cV^ {ci + Cr) 



Al^R 
A 

a: 



Cl + Cr cl'^ [Cl + Cr) 

Cl - Cr iCi {Ci - Cr) Zi Cl {ci - Cr) {2cf - 3cf Cr + 2cicl + cf) Z^ 



Cl + Cr c2 (Ci + Cr) 

I Cl [ci - Cr) Zi {Cl - Cr) zf 



(1.62) 

2 

^1.63) 



£'!_V^(C/ +C,,) 2DLcr^^ {Cl + Cr) 

Cl {-Cl + Cr) Zi ^ cf {Cl - Cr) zf 
Dlcl'^{ci+Cr) 2DlcV^ {Cl ^ Cr) 



Acf. {Cl + Cr) 

[-c,2+z(q2+c,2-qc,)] (1.64) 
[ci+i{2cr-ci)] . (1.65) 



In analogy to what was done for the "in" basis, this "out" basis can be used 
to obtain a decomposition of the (j) field operator as 
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did 



^v.outiout/i \ , --u, out tout / 1 \ I -v.outi out^/i _\ I -u,ouit out^ fi \ 

(1.66) 



in terms of the bosonic annihilation and creation operators for the out-going 
modes. This also leads to an alternative vacuum state defined by the condi- 



tions a'" 



Hilbert space. 



\0,out) 



^|0, out) = and an alternative "out" basis of the 



1.4-3 Bogoliubov transformation 

As both the "in" and the "out" basis are complete, the "in" and "out" scat- 
tering modes can be related by the simple linear scattering relations 

C"=^'C/+r'C*' (1-67) 
that can be summarized in terms of a unitary 2x2 scattering matrix 5* 

(1.68) 

Analogous relations hold for the (p^^ modes. 

Expressed in terms of mode amplitudes, these scattering relations define 
a linear Bogoliubov transformation relating the annihilation and creation 
operators for the "out" modes to the ones of the "in" modes. In the specific 
case of the sub-sub interface considered in the present section, all u and v 
modes involved in the scattering process have positive norm, so there is no 
mixing of the annihilation and creation operators: 

al^°^^ = Ral''" + T'ay. (1.69) 

As a result, the Bogoliubov transformation trivially reduces to a unitary 
transformation of the "in" and "out" Hilbert space that conserves the number 
of excitations and, in particular, preserves the vacuum state: if the system is 
initially in the |0, in) state with no incoming particles, the number of outgoing 
particles will also be zero, 

^v{u)^out ^ (o,m|a^(")'°"*^a^(")^™*|0,Mi) = 

= {Q,in\{T*{R*)al]'"'^+R'\T'*)al'"''^)^^ = 0. 

(1.70) 



No phonon can be created, but only scattered at the horizon. 
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Correlation functions are a modern powerful tool to investigate the properties 
of strongly correlated atomic gases [32] ■ We shall concentrate our attention 
to the correlation pattern of the density fluctuations at equal time, defined 
as 

G^^'>it;x,x') = lim(in|{ni(t,x),fii(i',a;')}|in) , (1.71) 

where {, } denotes the anticommutator. In our configuration with a spatially 
uniform condensate density n, the density fluctuation operator h\ can be 
expanded in the in-going annihilation and creation operators as 

= n{(f){x,t) + (t)\x,t)) = 

= n duj + Vl'^r) + + <i) + h-c] . (1.72) 

"'0 

or, alternatively, one can use the "out" basis. 

In either case, by evaluating G'^-' on the vacuum state |0, m) = \0,out), 
one flnds for one point located to the left (x < 0) and the other to the right 
{x > 0) the not too signiflcative expression 



G'^^\t;x,x') 



2Trmn{cr + q) 



{vo + ci){vo 



vo+ci 



va+Cr 



(1.73) 



which just shows correlations decreasing with the square of the distance 
weighted by the effective speed of sound in the different regions. The phys- 
ical origin of these correlations is traced back to the repulsive interactions 
between particles in the gas. 



1.5 Subsonic-supersonic configuration 
1.5.1 the modes and the matching matrix 

The warm up exercise discussed in detail in the previous section has allowed 
to take confidence with the formalism. In this section, we shall consider the 
much more interesting case of the acoustic black hole configuration sketched 
in Fig ll.3l taking again the fiow to be in the negative x direction {vq < 0), 
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Fig. 1.3 Upper panel: sketch of the subsonic-supersonic flow configuration. Low 
panels: dispersion relation of Bogoliubov modes in the asymptotic regions away from 
the horizon. 



we assume that the flow is sub-sonic Cr > \vo\ in the upstream a; > sector, 
while it is super-sonic c/ < \vo\ in the downstream a; < sector . 

The analogy with a gravitational black hole is simply understood: long 
wavelength sound waves in the a; < supersonic region are dragged away by 
the flow and no longer able to propagate in the upstream direction. The outer 
boundary of the super-sonic region separating it from the sub-sonic one plays 
the role of the horizon: long wavelength sound waves can cross it only in the 
direction of the flow, and eventually get trapped inside the acoustic black 
hole . Even if this picture perfectly captures the dynamics of long wavelength 
Bogoliubov waves in the sonic window where the dispersion has the form 
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Eq. (|1.32p . the supersonic correction that is visible in Eq. (|1.3ip introduces 
remarkable new effetcs as we shall see in the following of the section. 

The analysis of the dispersion relation and of the modes in the x > 
subsonic region on the right of the horizon is the same as given in the previous 
section: two oscillating modes exist with real wave vectors , s as well as 
two complex conjugate evanescent modes with fc!j-. 

In the a; < supersonic region on the left of the horizon, the dispersion 
relation has a significantly different shape, as shown in the lower-left panel 
Fig |1.3l In particular, it is immediate to see that there exists a threshold 
frequency LLimax above which the situation resembles the one of the sub- 
sonic regime: two oscillatory modes exist propagating in the downstream 
and upstream directions, respectively. Note that the upstream propagation 
occurs in spite of the super-sonic character of the underlying flow because 
of the super-luminal dispersion of Bogoliubov waves predicted by Eq. (|1.3ip . 
Of course, this mode falls well outside the sonic region where the hydrody- 
namic approximation is valid. The threshold frequency uJmax is given by the 
maximum frequency of the negative norm Bogoliubov mode as indicated in 
the lower-left panel Fig ll.3l In formulas, it corresponds to the Bogoliubov 
frequency of the mode at a kmax value such that 



1 

Ti 




1/2 



(1.74) 



The < Lu < LOmax casc is much more interesting: from Fig. II. 3[ one sees 
that four real roots of the dispersion relation exist, corresponding to four 
oscillatory modes, two on the positive norm branch and two on the negative 
norm one. Two of these modes denoted as u, v lie in the small k region at 



kv 
kii. 



V - Ci 

V + Cl 



-I ''I 



8(«o - cif 

cfzf 
Sivo + cif 



+ Oizf) 



+ Oizf) 



(1.75) 
(1.76) 



and have a hydrodynamic character. Differently from the sub-sonic case, both 
of them propagate in the downstream direction with a negative group velocity 
^ < 0: also the u mode that in the comoving frame with the fluid propa- 
gates to the right is dragged by the super-sonic flow and turns out to be 
forced to propagate in the left direction. Furthermore, while the k^, solution 
belongs as before to the positive norm branch, the fc„ solution belongs now 
to the negative norm branch and the corresponding excitation quanta carry 
a negative energy cu < 0. 

The wavevector of the other two roots indicated as k^ and k^ in the figure 
is non-perturbative in ^ 
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ujvo 



0{zf 



± 



ciii 



icf + 2vi)cfzt 



+ 0{zf) 



(1.77) 



and lies well outside the hydrodynamic region. Comparing these roots with 
Eq. (|1.39p . one realizes that fcs^ arc the analytic continuation for supersonic 
flow of the growing and decaying modes previously discussed for the subsonic 
regime. The mode belongs to the positive norm branch, while to the 
negative one; both of them have a positive group velocity and propagate in 
the upstream direction. 

For u} < LUmax, the general solution of the modes equation in the super- 
sonic (left) region reads then 



I rpl Ak^,x 



while in the sub-sonic (right) region it reads 



A'',D\e'''+''+ALDLe'' 



(1.78) 

An analogous expression holds for ipl^ once we replace D[lij) with E{ll>). 

As we have discussed in the sub-sub case, the field amplitudes on the left 
and the right of the discontinuity point at x = are related by 



A'- 

4 



M 



\ALJ 



(1.79) 



the matching matrix being written as M — Wr in terms of Wr given by 
Eq. prig)) and 



Wi 



( D 

i 

\ikiEi iklEl ik^Dl ik{D\J 



iklDl iklDl ik'^D\ ik{D\ 



U U 



(1.80) 



1.5.2 The "in" and "out" basis 



We can now proceed to construct the "in" scattering basis. Differently from 
the sub-sub case discussed in the previous section, there are now three "in" 
scattering modes associated to the processes sketched in Fig |1.4l 



1 Simple models of Hawking radiation in atomic BECs 



25 






Fig. 1.4 Sketch of the Bogoliubov modes involved in the "in" (left panels) and "out" 
(right panels) basis. The mode labels refer to the dispersion shown in the lower panels 
of FigO 



The mode (/>Ji''" is defined as an initial left-moving unit amplitude v wave 
propagating in the x > sub-sonic region towards the horizon, which upon 
scattering generates in the subsonic region a reflected right-moving u wave 
of amplitude and a spatially decaying wave of amplitude . The trans- 
mitted waves in the x < supersonic region arc now in the number of two, 
and both travel in the leftward direction along the flow. One is the standard 
transmitted v wave, with positive norm and amplitude the other is a neg- 
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ativc norm u wave with amplitude A^, the so-cahed anomalous transmitted 
wave. 

To leading order in zi, the corresponding amplitudes are 



Cr Va - Cl 



Cl Vq ~ C 
VQ + Cr _ 
Vq - Cr 



S'U: 



Lvr: 



Cl Cr 



Lvr: 



Cl-J^l^J Crivl - C\) 



V2i?;(i>o-Q)(c2-v2)3/2(c, + cO 



(1.81) 
(1.82) 
(1.83) 

y 

(1.84) 



Note the shorthand notation introduced in (|1.8mi.8^ to simply identify the 
incoming and outgoing channel: for example the matrix element Sui,vr indi- 
cates that the incoming channel (second index) is a w-mode entering from the 
right region, while the outgoing channel (first index) is a u-mode escaping 
in the left region. The conservation of the Bogoliubov norm translates into a 
unitary condition between the amplitudes of the propagating modes, 



\A. 



I |2 



\A' 



I |2 



1, 



(1.85) 



where the minus sign comes from the negative norm u, Z-mode. 

As it is sketched in Fig ll.4[ the other two "in" scattering modes (l>l]™ and 
^4,in constructed in a similar way by imposing a unit amplitude in the 
/c3 or k4 waves incident on the horizon from the left supersonic side. For the 
03,2n itjj^!> scattering mode, the corresponding amplitudes are given by 



{vi-cf)yHvo+' 

^3/2 



Al 



c/ V2zi{ci + Cr)\Jcl - v\ 
V2^(«g-C;2)3/4(^;o+e,) 
Q^(c2 - c\)^ cl - vl 
(^Vl-C}fl\v^+Cr) 
C^^ ',/2zi{ci - Cr)^cl - vl 

:2— 2Ti^'" 



2 - 



• 2Dl{vt,-ct., 
and satisfy the unitarity relation 



cl-vl) ^ S+r,3l. 



Aif = l. 



<5'«r,3/(l-87) 

■ •S'ui, 3(1-88) 
(1.89) 

(1.90) 
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As a most remarkable point, note that the amphtudes for the propagating 
modes now diverge at small lo as -k=. 

For the "in" scattering mode one gets instead 

(»ii-cf)=/-'(»„+c,) 



/2^{vi~cff/'^{va + Cr) 



and the unitarity condition reads 

\Al\'' + \Al\^-\Ai\^ = -l: (1.95) 

the minus sign on the right-hand side comes from the fact that the incoming 
unit amplitude fc4 mode has negative norm. All together, the w, 3 and 4 "in" 
scattering modes form a basis on which to expand the </> field operator. 

The construction of the "out" basis proceeds along similar lines: one has 
the three i/)""*, (t)'^^ ^^i^ "out" scattering modes, where the l{r) label 
near the superscript u indicates again the left (right) region of space. The 
corresponding scattering processes are depicted in Fig ll.4l 

As the corresponding amplitudes will not ne needed in the following, we 
refer the reader to Ref. [T7| for their explicit expression. As discussed in [6l [18] , 
the field operator can then be equivalently expanded in the basis of the "in" 
scattering modes as 

or, equivalently on the basis of the "out" scattering ones. Note in particular 
the third term on the right-hand side, as the corresponding mode is a 
negative norm one, this term enters with a c/)™; field multiplied by a creation 
^4,mt operator: as we shall see in the next sub-section, this simple fact is 
the key element leading to the emission of analog Hawking radiation by the 
horizon. 
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1.5.3 Bogoliubov transformation 

The "in" and "out" basis are now related by a 3 x 3 scattering matrix 5* 
relating the three incoming states to the three outgoing states. Explicitly 



(1.97) 

'?^3" = Svl,3l(f>v'!'l + Sur,3l4'Tr + Sul,3l(t>u"l : (1.98) 
"^4" = Svl,4l(j)y'^i + SurAl4'u'!r + SulAl'P'^l ■ (1.99) 

Because of the negative norm of the (/)°"* mode, conservation of the Bogoli- 
ubov norm imposes the modified unitarity condition S'^rj S = Srj with 



77 = diag{l, 1,-1) and the scattering matrix S mixes positive and negative 
norm modes. 

As a result, the Bogoliubov transformation relating the creation and de- 
struction operators of the "in" and "out" scattering states is no longer trivial 
and mixes creation and destruction operators as follows 

Svl.vr Syl_3l SylAl \ I a^'™ \ 

Sur,vr Sur,3l SurM . (1.100) 

SuLvr SuLZl Sul,Al / \ O^*" V 

The non triviality of the Bogoliubov transformation has the crucial conse- 
quence that the "in" and "out" vacua no longer coincide |0, m) ^ |0,owt): 
while the |0, in) "in" vacuum state (defined as the state annihilated by 
the a^j"'^'*^''" operators) contains no incoming phonons, it contains a finite 
amount of phonons in all three out-going modes due to a parametric conver- 
sion process taking place at the horizon, 

= (0,*n|ar°"*^ar°"*|0,«ri) = \SurM? , (1-101) 
nl^' = (0,*n|a:i''°"*ta^'^-*|0,zn) = + \SuL31? . (1-102) 

<f = (0,*7i|ai;^™*ta-,-*|o,zn) = |5„i.4ip. (1-103) 

Note in particular the remarkable relation 

K'' = \Sul,vr\'' + \Sul,3l? = \SurAl? + - + n"J . (1.104) 

The physical meaning of the above relations can be understood as follows. 
Suppose that at t = — oo we have prepared the system in the |0,m) vacuum 
state, so there are no incoming phonons. We are working in the Heisenberg 
picture of Quantum Mechanics, so that |0,m) describes the state of our sys- 
tems at all time. Now eqs. (|1.101j|1.103|) tell us that at late time in this state 
there will be outgoing quanta on both sides of the horizon: the vacuum has 
spontaneously emitted phonons. This occurs by converting vacuum fluctu- 
ation of the fc4 mode into real on shell Bogoliubov phonons (sec lower left 
panel of Fig ll.4p in the hydrodynamic region. 
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While processes involving particle creation in time-varying settings are 
well-known in quantum mechanics, e.g. the dynamical Casimir effect P^[2D] . 
the production of particles in a stationary background seems to contradict 
energy conservation . The solution of this puzzle relies in Eq. (jl.l04l) : besides 
the positive energy ur and v phonons , there is also production of negative 
energy (ul) phonons, the so called "partners" which propagate down in the 
supersonic region. The number of these latter equals the number of the for- 
mers. This is how energy conservation and particles production coexists in 
our stationary systems. As particles are produced in pairs with opposite ±io 
frequencies, energy is conserved. 

Now let us give a closer look at Eq. (|1.101[) : according to this, an hypo- 
thetical observer sitting far away from the horizon in the subsonic region at 
X — > +00 will reveal a flux of phonons coming from the horizon. This is just 
the analogue of Hawking black hole radiation . The number of phonons of 
this kind emitted per unit time and per unit bandwidth is 

The ^ behavior of the above expression is reminiscent of the low frequency 
expansion of a thermal Bose distribution jl8j 

- 1 

and one can try to identify the — coefficient of p.lOSp as an effective tem- 
perature 

^ _ h iCr+ Vq) {v^ - Cf 2Cr ^^^^ 

ks [Cr - Vo) {cl - cf) Q^i 

As the surface gravity of our toy model with an abrupt discontinuity in the 
flow is formally infinite while the temperature remains finite, the connection 
of the analog model to the original gravitational framework seems to fail. 
However, to investigate the correspondence with the gravitational black holes, 
one has to consider more general and realistic velocity profiles where the 
transition from the subsonic region to the supersonic one is smooth enough to 
justify the hydrodynamical approximation. Accurate numerical calculations 
in this regime show that the emission is indeed thermal in this case and the 
temperature is to a good accuracy determined by the surface gravity k of the 
associated black hole according to Eq. (|1.2ip . As it was shown in [SJ , the 
original Hawking's prediction for the emission temperature holds provided 
the spatial variation of the flow parameters occurs on a characteristic length 
scale longer than ^^/'^k"^/'^. Of course, the thermal spectrum is restricted 
to frequencies lower than the upper cut-off at uJmax- above this frequency, 
one in fact recovers the physics of the sub-sub interface where no emission 
takes place. These results, together with the full numerical simulation of [H] 
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confirm that the emission of Hawking radiation is not an artifact of the 
hydrodynamical approximation and provide an independent vahdation of the 
model of phonon propagation based on the metric Eq. (|1.20p . 

Even in the most favorable configurations, realistic estimates of the Hawk- 
ing temperature in atomic BECs give values of the order of 10 nK, that is 
one order of magnitude lower than the typical temperature of the condensates 
(100 nK). This makes the Hawking emission of Bogoliubov phonons in BECs 
a quite difficult effect to reveal in an actual experiment, as the interesting 
signal is masked by an overwhelming thermal noise. 

A proposal to overcome this difficulty was put forward in [11]: as the 
pairs of Bogoliubov excitations produced by the Hawking process originate 
from the same vacuum fluctuation, their strong correlation is expected to be 
responsible for specific features in the correlation function of density fluctua- 
tions. This idea was soon confirmed by numerical simulations of the dynamics 
of atomic condensates in acoustic black hole configurations. The features ana- 
lytically predicted in |llj are indeed visible in the density correlation pattern 
and, moreover, are robust with respect to a finite temperature. At present, 
this method represents the most promising strategy to experimentally detect 
the analog Hawking effect in atomic Bose-Einstein condensates . A first in- 
vestigation of the power of density correlation techniques in the context of an 
analog dynamical Casimir effect in condensates has been recently reported 
in [22] along the lines of the theoretical proposal in [21] . 



1.5.4 Density- density correlations 

The fact that the density correlation function in BECs exhibits characteristic 
peaks associated to the phonons creation a la Hawking can be easily seen in 
our simple toy model in the following way. For simplicity, let us restrict our 
attention to the contribution to the density-density correlation function due 
to the "out" particles and consider the decomposition 

-l/i _\ ^ , _ / ["V.out/iout I out\ I I ^ur,out/iout , out\ , 

Jo 

+ aZ''''-'\ruJ+vZ') + h.c.]. (1.108) 

Expanding the "out" creation and annihilation operators in terms of the "in" 
ones and using the approximate form of the S matrix elements given by Eq. 
(|1.100p . and finally evaluating expectation values on the |0, in) "in" vacuum 
state, one finds that the above expression describes correlations between the 
(ur) and (ul) particles and the (ur) and (vl) particles if the points x and 
x' are taken on opposite sides with respect the horizon, while one finds {ul)- 
(vl) correlations if both points are inside the horizon. If both x, x' are located 
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Fig. 1.5 Color plots of the rescaled density correlation function (nofr) x 
[G^'^\x,x') — 1] a time grUt/h = 160 after the switch-on of the black hole hori- 
zon. The calculation has been performed using the truncated- Wigner method of | 12 |. 
Black hole parameters: |i!o|/ci = 1.5, luol/cr = 0.75. The dashed white line indi- 
cates the analytically expected position (|1.110|l of the negative peak in the density 
correlation signal. 



outside the horizon, correlations just show a monotonia decrease with distance 
as in the sub-sub case. 

As in general one has |S't,i.4i| <C |S'„;.4;|; the main contribution to the 
density correlation in the a; < and x' < sector comes from the {ul)- 
(ur) term describing correlations between the Hawking phonon (ur) and its 
partner (ul). Integrating over all frequencies upto Umax, one obtains term of 
the form [111 



1 



r.2\3/2 



sm 



+ Ci > 



vo+ci 



Aim Cl{vo + Ci){vo - Cr){Cr - Ci) 



VO + Cr 



vo+ci 



(1.109) 

From this expression, it is easy to see that the density-density correlation 
function has a negative value and is peaked along the half-line 
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^— = . (1.110) 

Vo + Cr Vq + Cl 

The stationarity of the Hawking process is apparent in the fact that the peak 
value of Ea. (|1.109p does not depend on the distance from the horizon. 

The physical picture that emerges from this mathematical derivation is 
that pairs of (ul) and (ur) phonons are continuously created by the horizon 
at each time t and then propagate on opposite directions at speeds Vui = 
Vo + Cl < for (ul) and Vur = vq + Cr > for (ur). At time At after their 
emission they are located at x = Vui At and x' = Vur At, which explains 
the geometrical shape of the peak line Eq. (jl.llOp where correlations are 
strongest. An example of numerical picture of the correlation function of 
density fluctuations is shown in Fig ll.51 the dashed line indicate the expected 
position of the peak line Eq. (|1.110p . A detailed discussion of the other peaks 
(that are barely visible on the color scale of the figure) can be found in p^lTS] . 



1.5.5 Remarks 

Let us try to summarize the results discussed in the present chapter. We 
have seen that if a stationary flowing BEC shows a horizon-like boundary 
separating an upstream subsonic region from a downstream supersonic one, 
a spontaneous emission of Bogoliubov phonons occurs at the horizon by con- 
verting zero-point quantum fluctuations into real and observable radiation 
quanta. The emitted radiation appears to an observer outside the horizon in 
the subsonic region to have an approximately thermal distribution: this is 
the analog Hawking effect in BECs. 

As in a gravitational context nothing can travel faster than light, the hori- 
zon has a well-defined meaning of surface of no return: no physical signal 
can travel from inside the black hole to the outside crossing the horizon in 
the outward direction. In the case of an atomic BECs, the "sonic" horizon is 
defined as the surface where the speed of sound c equals the velocity |uo| of 
the fluid: for what concerns the hydrodynamical w, v modes at low wavevec- 
tor, the acoustic black hole exactly mimics what happens in gravity: no long 
wavevector sound wave can cross the horizon in the upstream direction. On 
the other hand the dispersion of Bogoliubov modes in an atomic BEC shows 
significant super-luminal corrections: the higher the wavevector of the exci- 
tation, the larger its group velocity. As a result, the /c3,4 modes are able to 
travel in the upstream direction in the super-sonic region inside the horizon 
and therefore to escape from the black hole. In contrast to the hydrodynamic 
modes, they are not trapped inside and do not see any horizon: for them the 
gravitational analogy has no meaning. 

Some authors have recently introduced wavevector dependent rainbow 
metrics to describe the propagation of different modes at different wave vec- 
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tors and have defined several distinct concepts of horizon, such as the phase 
horizon and the group horizon. Our opinion is that these additional concepts 
may end up hiding the essence of Hawking radiation behind unessential de- 
tails. 

The key ingredient in order to have the emission of radiation in the "in" 
vacuum state is in fact the presence of negative energy states that allow to 
emit a pair of quanta while conserving energy: this requires that the flow 
undergoes supersonic motion in some spatial region. The main role of a hori- 
zon where the flow goes from sub- to super-sonic is to determine the thermal 
shape spectral distribution of the emitted radiation. To better appreciate this 
fundamental point, the next chapter will be devoted to a short discussion of 
configurations with a super-sonic flow on both sides of the interface: as it was 
first pointed out in |26| . a spontaneous emission of radiation takes place in 
this case in spite of the total absence of a horizon: sound waves are always 
dragged by the super-sonic flow and can not propagate upstream. Because of 
the absence of an horizon, the resulting spectral distribution of the associ- 
ated zero-point radiation is however very different from the thermal Hawking 
radiation, with a low-frequency tail dominated by a constant term instead of 



1.6 Supersonic-supersonic configuration 

Consider a BEC undergoing an everywhere supersonic motion , with a sound 
velocity profile varying abruptly at a; = 0: no sonic horizon is present in this 
setting and at all points long wavelength sound waves are dragged in the 
downstream direction by the underlying flowing fluid. 

The dispersion relation pattern on either sides of the discontinuity is shown 
in the lower panels of Fig |1.6l For u < uj„iax = inin[w^^^, w^^^,], one has four 
oscillatory solutions in both regions with real wavevectors. The fc„ and fc„ 
hydrodynamic solutions propagate in the downstream direction (i.e. to the 
left with negative vo ) while the large wavevector ks and ki solutions are able 
to propagate upstream. While the fc«,3 solutions correspond to positive norm 
modes, the ku,i are negative norm ones. 

The general solution of the mode equations in both regions reads 




'Mr) r)l{r) rk'^x , Mr) Tjl{r) tk[ 



X 



— e 



+ 



(1.111) 



As usual, the left and right amplitudes are related by 
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Fig. 1.6 Upper panel: sketch of the supersonic-supersonic flow configuration. Low 
panels: dispersion relation of Bogoliubov modes in the asymptotic regions away from 
the horizon. 



M 



the matching matrix being given by M = W'^Wr with 



/ 



r,i{r) 



D 



l{r) 



D 



l{r) 
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I(r) 



(1.112) 



(1.113) 
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As sketched in Fig ll.7[ the "in" basis is here defined by four in-going waves: 
two of them {u, v) are incident on the discontinuity from the right (left panels 
on the first and second rows); the two others (3,4) are incoming from the left 
(left panels on the third and fourth rows). The "out" basis is defined along 
the same lines as sketched in the four panels of the right column. 

The field operator can be expanded either in the "in" or in the "out" basis 

as 

(/)=/ duj[aj (p^^r + aj '(Pu.r+aj 93.1 + aj '04,;+ 
Jo 

+ aj Vi,,r+aj '^u,r+aj '^3,1 +au: ^4,1] (1-114) 

or 

Jo 

+ aj +aj 'Pu,i + « J V3.r + a J f4,r\- (1-115) 

The relation between the "in" and "out" basis are 

— Svl^vr4>v^l + Sul,vr<j>'^l + S^r.vr<t>'3^r + '5'4r,i;r04"* , (1.116) 

= Svl^3l<^°"l + SuLSl't'Z^l + S3r,3l(f>3"r + S4r,3l(f>4"r i (1.118) 
^4^1 = '5'i>Z,4/0°"* + Sul,4l<i>°i!^l + SzrAl^S^r + 5'4r,4!04"* (1.119) 

and the corresponding relation between the "in" and "out" annihilation and 
creation operators reads 









Svl,ur 


Svl,3l 


Svl,4l \ 




1 


^ \ 








Sulfur 


Sul,3l 


Sul,4l 








p.3,out 




S3r,vr 


S3r,ur 


S3r,3l 


S3r,4l 












^ S4r,vr 


S4r,ur 


S4r,3l 


S4r,4l ) 




V 





(1.120) 



Explicit expressions for the corresponding 16 amplitudes are listed in the 
Appendix. 

Here we see again that the S matrix mixes creation and annihilation op- 
erators. As a consequence, the |0,m) and |0,oui) vacua do not coincide: in 
particular the "in" vacuum |0,m) state with no incident quanta leads to a 
finite amount of out-going particles that can be detected: in the left region, 
they belong to the u, v modes; in the right region, they belong to the 3, 4 
modes. 

More precisely 
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Fig. 1.7 Sketch of the Bogoliubov modes involved in the "in" (left panels) and "out" 
(right panels) basis. The mode labels refer to the dispersion shown in the lower panels 
of FiglTel 
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n 



v,l _ 



(0,zn|a3'°"*ta3. 
(0,*n|ai'-*ta4. 




|0,m) = \S,l,ur\^ + \S,lAl\^ 

|0,zn) = \S3r,ur\^ + \S3rAl\^ 
\0,in) = \S4r,vr\'^ + 154^,3;^ 



(1.121) 
(1.122) 
(1.123) 
(1.124) 



'UJ 



u,l 



n] 



n, 



and unitarity of the S matrix imposes that 




+ n. 



4,r . 



(1.125) 



the number of positive energy particles equals the number of negative energy 
ones. From the explicit expressions for the S matrix elements listed in the 
Appendix, it is immediate to see that all spectral distributions n^'', n"'', 
nj^''' and nfj''" at low frequencies are dominated by constant terms, in stark 
contrast with the 1 /w shape of the thermal Hawking radiation. 



1.7 Conclusions 

In this chapter, we have given an introductory review to Hawking radiation 
effects in atomic Bose-Einstein condensates Bose-Einstein condensates. Fo- 
cussing our attention on a simple toy model based on a piecewise uniform 
flow interrupted by sharp interfaces, we have made use of the standard Bo- 
goliubov theory Bogoliubov theory of dilute condensates to obtain analytical 
predictions for the quantum vacuum emission of phonons that is emitted by 
the interface: necessary and sufficient condition for this emission to occur 
is that the flow be somewhere super-sonic. While the low- frequency part of 
the emission follows an approximately thermal form for a black-hole inter- 
face separating a sub-sonic upstream region from a super-sonic downstream 
one. a completely different spectrum is found for flows that do not show any 
horizon and are everywhere super-sonic. 

The interest of our development is manyfold: on one hand, our analytical 
treatment provides an intuitive understanding of Hawking radiation based on 
a Bogoliubov generalization of the scattering of waves by square potentials in 
one dimensional Schrodinger equation. On the other hand, our derivation is 
however completely "ab initio" , based on the fundamental microscopic quan- 
tum description of the BEG without any recourse to the gravitational analogy. 
As a result, it does not depend on the hydrodynamic approximation that un- 
derlies the introduction of the effective metric and shows that the Hawking 
effect in atomic BECs is not at all an artifact of the low wavelength (hydro- 
dynamical) approximation: no transplanckian problem is present which may 
cast doubts on the derivation, rather our derivation shows that the trans- 
planckian problem is itself an artifact of the hydrodynamical approximation. 

The intense theoretical and experimental activity that is currently in 
progress makes us confident that the existence of analog Hawking radiation 
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will be soon experimentally confirmed. The robustness of Hawking radia- 
tion with respect to the microscopic details of the condensed-matter system 
would be a strong indication that, in spite we do not have any knowledge 
of the quantum microscopic description of gravity, Hawking's prediction of 
black hole radiation with its important thermodynamical implications is a 
real milestone in our understanding of Nature. 

Acknowledgements Continuous stimulating discussions with S. Finazzi, R. Parentani 
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from ERG via the QGBE grant. 



Appendix 

In this Appendix, we give the explicit expressions for the leading order in 
the small w limit of the S Matrix coefficients for the supersonic-supersonic 
configuration treated in Sect. 11.61 Note in particular how the ones involved 
in the vacuum emission {Svi^ur, Sui.vr, S^r.Ah 5'4r,3/) grow as ^/uJ at low lu, 
while {S4r,vr, S3r,ur, SuL3i, Syi^n) tend to constant. 
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